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I. INTRODUCTION 

Relaxation dynamics of random field Ising model (RFIM) at zero temper- 
ature provides a simple caricature of complex non-equilibrium phenomena. 
RFIM with ferromagnetic as well as anti-ferromagnetic interactions has been 
investigated in the context of two distinct classes of relaxational behavior. 
The ferromagnetic RFIM shows relaxation by avalanches and has been ap- 
plied []I],|2| to study hysteresis and Barkhausen jumps in magnetic materials. 
The anti-ferromagnetic model does not support avalanches, and is more ap- 
propriate [^,^ for the glassy kind of dynamics where relaxation may proceed 
by single localized events. In spite of the simplicity of RFIM and its relax- 
ational dynamics at zero temperature, exact solution of the model is difficult 
on account of its randomness. So far, analytic solution of the model has been 
obtained in one dimension [j5[] and on a Bethe lattice [^,0 for ferromagnetic 
interactions only. In this paper, we present a solution of the non-equilibrium 
dynamics of the anti-ferromagnetic RFIM at zero temperature in one dimen- 
sion in a limited range of the applied field. For simplicity, we limit ourselves 
to a bounded uniform distribution of the random field. It is possible to ex- 
tend this solution to other distributions of the random field, but this will not 
be taken up in this paper. 



II. THE MODEL 

At each site of a one dimensional lattice, there is an Ising spin Si = ±1, 
i = 1,2,3 , n which interacts with its nearest neighbors through an anti- 
ferromagnetic interaction J. Each site has a quenched random field hi drawn 
independently from a continuous bounded distribution, 

p{h,) = 2^ if -A<h,<A 

= 0, otherwise. (1) 

The entire system is placed in an externally applied uniform field ha- The 
Hamiltonian of the system is given by 

H = -JJ2stS,+i-J2htS^-haJ2st (2) 

i i 

We consider the hysteretic response of this system when the external field 
ha is slowly increased from — oo to +oo. We assume the dynamics to be the 
single-spin-flip Glauber dynamics at zero temperature, i.e. a spin is flipped 
only if it lowers the energy. We assume that if the spin-flip is allowed, it occurs 
with a rate F, which is much larger than the rate at which the magnetic field 
ha is increased. Thus we assume that all flipable spins relax instantly, and 
the spin Si always has the same sign as the net local field k at the site. 

Sj = sign U = sign [J(sj_i + s^+i) +h^ + ha] (3) 

The hysteretic response of the system to an applied field ha is measured 
by the magnetization m{ha) per spin. 



m{ha) = -Y.s^ (4) 

We start with ha = — oo, when m = — 1, and increase ha slowly to ha = +00 
when m=l. Our object is to calculate m{ha) for all values of ha- The 
magnetization mR{ha) on the return trajectory when ha is slowly decreased 
from ha = +00 to ha = —00 can be obtained from m{ha) by a symmetry 
relation mR{ha) = —m{—ha). Therefore the knowledge of the magnetization 
curve on the lower half of the hysteresis loop suffices to determine the entire 
hysteresis loop, and we can limit ourselves to the calculation of m{ha) alone. 

III. SIMULATIONS 

It is instructive to look at a computer simulation of the model before 
proceeding to obtain it analytically. Figure 1 shows magnetization m{ha) 
in an increasing applied field ha for an anti-ferromagnetic RFIM (J = —1), 
obtained from a simulation of 10^ spins averaged over 10^ different realiza- 
tions of the random field distribution of width 2 A = 1. The magnetization 
m{ha) rises from -1 to +1 in three steps. We call these steps as ramp-I 
{ha = -2\J\ - A to ha = -2\J\ + A); ramp-II {ha = -A to ha = +A); 
and ramp-III {ha = 2| J| — A to /ia = 2| J| + A). The ramps are connected 
to each other by two plateaus; plateau-I { ha = — 2| J| + A to /ia = —A); 
and plateau-II {ha = +A to ha = 2\J\ — A). On the plateaus, the magneti- 
zation remains constant even though the applied field continues to increase. 
Plateaus occur for A < |J| (small disorder), and simulations suggest that 



magnetization on the plateaus is independent of A. Numerically, the mag- 
netization on the plateaus is approximately m^ = —.135 on plateau-I, and 
m^^ = .109 on plateau-II. 

The qualitative shape of m{ha) is easy to understand. Due to the anti- 
ferromagnetic interaction between nearest neighbors, spins with both neigh- 
bors down are the easiest to be turned up by an applied field increasing in 
the up direction. Such spins turn up on ramp-I. Next are the spins with 
one neighbor up and one down which turn up on ramp-II. Spins with both 
neighbors up require the largest applied field to turn up, and these turn up 
on ramp-III. On each ramp, the sequence in which the spins turn up is de- 
termined by the distribution of the quenched random field. Spins with large 
positive quenched field turn up before spins with a lower quenched field. The 
quenched field lies in the range —A to +A. Thus each ramp has a width 2A 
along the axis of the applied field. 

When a spin turns up on ramp-I, its nearest neighbors are placed in a 
category so that they cannot turn up before ramp-II. Similarly when a spin 
turns up on ramp-II, its nearest neighbor which is down cannot turn up 
before ramp-III. This is essentially the reason for the absence of avalanches 
in the anti-ferromagnetic RFIM. Occasionally on ramp-II and ramp-III, a 
spin turning up can induce its nearest neighbor which is already up to turn 
down. We will discuss the situation on ramp-II in detail later, but suffice 
it to say here that this process too does not cause an avalanche. With 



anti-ferromagnetic interactions, spins turn up one at a time (no Barkhausen 
noise), and the calculation of m{ha) becomes essentially a matter of sorting 
quenched random fields in decreasing order on each ramp. The difficulty 
arises from the fact that the a posteriori distribution of random fields on 
spins classified according to the ramp on which they turn up is significantly 
modified from the uniform distribution given in Equation (1). Indeed, the 
main object of this paper is to calculate this distribution. 

IV. RAMP-I 

Magnetization on ramp-I was determined earlier in reference 5 by exploit- 
ing a similarity between this problem and the problem of random sequential 
adsorption(RSA) [[Tj. The rate equations of the RSA problem were used to 
determine m{ha) on ramp-I, but they could not determine m{ha) on ramp-II 
and ramp-III. Here we rederive the result for ramp-I by a different approach 
which can be extended to ramp-II as well. 

We start with a large negative applied field {ha = — cxd) when all spins are 
down and increase the applied field slowly. At the start of ramp-I, i.e. at 
ha = — 2| J| — A, the spin with the largest positive quenched field becomes 
unstable and flips up. As the appHed field continues to increase, the spin 
with the next largest quenched field turns up unless it happens to be next 
to a spin which is already up. In this case the spin which turns up next 
is the one with the largest quenched field from among the spins with both 



neighbors down. Consider an arbitrary point on ramp-I at an applied field 
ha = —2\J\ — h. In the following, the field h will be used more frequently than 
ha- In general, the field h will be given by the relation h = —ha mod |2J|, so 
that on each ramp it has the same range as the random field (—A < h < A). 
At an applied field ha = — 2|J| — h all spins with quenched random field 
hi > h are relaxed, and a fraction of them are up. The fraction of sites with 
hi > h is given by 

p{h) = j^ p{hi)dh = ^^ (5) 

The fraction of spins which are up on ramp-I at ha = — 2| J| — h is given 
by 

P/ = l[l-e-2^] (6) 

Let Pjj^ be the probability (per site) of finding a pair of adjacent down 
spins on ramp-I at the applied field ha = —2\J\ — h. It can be calculated as 
follows. Imagine coloring all sites with hi > h black, and all sites with hi < h 
white. Consider two adjacent down spins A and B shown in Figure (2). The 
sites A and B can be both white, both black, or mixed. Given that A is 
down, it is clear that the state of B can only be infiuenced by the evolution 
of the system to the right of B. Similarly, given that B is down, the state of 
A can only be infiuenced by the evolution of the system to the left of A. We 
shall refer to this as the principle of conditional independence [@] . It requires 

P/^ = P{A i \B i)P{B i \A i) (7) 

7 



where P{A j \B j) is the probabihty that spin at site A is down given that 
spin at B is down, and P{B I \A j) is the probabihty that B is down given 
that A is down. We take up the calculation of P{B I \A I). If B is a white 
site, P{B I \A I) = 1 because the white sites have not been relaxed from 
their initial state. If B is a black site and the site to the right of B is a white 
site then P{B I \A I) = 0. In general P{B I \A j) depends on the length of 
the string of black sites to the right of B. Suppose B is a black site, and there 
are (n-1) additional black sites to the right of B. In this case, the probability 
Pg that B is down satisfies the following recursion relation, 

Pb = -Pb-' + (1 - -)Pb-' (8) 

n n 

The rationale for the above recursion relation is as follows. Let the black site 
farthest from B on the right be labelled as the n-th site. Any of the n sites 
could flip first. The probability that the n-th site flips first is therefore equal 
to -. If this happens, (n-l)-th site is prevented from flipping up on ramp-I. 
The probability that B is down is now reduced to the probability that the 
end point of a chain of (n-2) black sites is down i.e. Pb ^. This accounts for 
the first term in equation (8). The probability that n-th site is not the first 
site to flip up is equal to (1 — -). Given this situation, the probability that 
B is down is equal to the probability that the end of a string of (n-1) black 
sites is down. This accounts for the second term in equation (8). We can 
rewrite the recursion relation (8) as 

{Pb - Pb-') = --[Pb-' - Pb-'] (9) 

n 



It has the solution, 

pe=j: ^—^ (10) 

Summing over various possible values of n with appropriate weight, we 
get 

oo n ( _'\\m 

n—Om—O ^• 

oo (_-\\m oo oo ( _^\m 

= T.^-^ii-p)T.P''= E^ = e-^ (11) 

Thus, 

P'u = ^"* (12) 

Let P|^ be the probability per site of finding a down spin and Pj^ the 
probability per site of finding a down spin which is followed by an up spin. 
Clearly, 

Pl = p'u + P'» = 1 - ^/ (13) 

Keeping in mind that on ramp-I an up spin must be preceeded (as well as 
followed) by a down spin, we get P|^| = P^ . Thus, 

P] = 1 - P|l - P| 



or. 



Pl = \[^-Pli] = \[^-^-''] (14) 



2' ^^'2 

The magnetization on ramp-I is given by 



m\h) = 2P/(/i) - 1 = -e-2M^) (15) 

Equation (15) has been superposed on the simulation data for ramp-I 
shown in Figure 1. The fit between the simulation and theory is so good 
that the two curves are indistinguishable from each other on the scale of 
Figure (1). The exact value of the magnetization on plateau-I is equal to 
— ^ which is approximately equal to -.135. 

V. PLATEAU-I 

Plateau-I contains down spins in singlets and doublets punctuated by up 
spins. The down spins were relaxed on ramp-I but did not turn up because 
the then applied field was not strong enough to turn up spins with at least 
one neighbor up. The singlets have two neighbors up, and the applied field 
on ramp-II is still not strong enough to turn them up. They have to wait for 
their turn on ramp-III. On the other hand, each spin in a doublet has one 
neighbor up and one down. It therefore experiences zero net field from its 
neighbors. The net field on a doublet spin is simply the sum of the random 
field hi on its site and the applied field ha = —h. It turns up when ha+hi > 0, 
or hi > h. Note that the random field is bound in the range —A < hi < +A. 
Therefore an applied field smaller than —A is sufficiently negative to pin 
down all doublets. This is the reason for the plateau in the magnetization 
for applied fields in the range (— 2|J| + A)</i<— A. In each doublet, the 
spin with the larger quenched field hi flips up on ramp-II when the applied 
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field reaches a value such that hi > h. The spin with the smaller quenched 
field then becomes a singlet which does not fiip up before ramp-III. Thus, in 
order to find the form of ramp-II, we need to find the a posteriori distribution 
of quenched random fields on the doublets. 

Consider a doublet on plateau-I as shown in Figure (3). The doublet sites 
are denoted as 1 and 2, and the quenched random fields on these sites are hi 
and /i2. The probability (per site) of finding a doublet on plateau-I is easily 
obtained from equation (12) by putting p = 1. It is equal to -^. We now 
calculate the probability distribution for hi and h2. The distribution of hi 
and /i2 will be identical if the evolution of the system on the two sides of 
the doublet is similar to each other. We assume this to be the case for now, 
although we shall examine a more general situation later. In order to obtain 
the desired probability distribution, consider the system on ramp-I when all 
spins are relaxed upto an arbitrary applied field —2\J\ — h. At this point, the 
probability per site for finding a doublet is Pi i = e~^^, where p is fraction of 
sites with hi > h (black sites) on an infinite lattice. 

It is often easier to think in terms of the fraction p. We therefore introduce 
similar fractions for the quenched fields given by the following relations. 



Pih^) = ^^; (i=l,2 



Given that site 1 is down on ramp-I, the conditional probability that the 
adjacent site 2 is also down is equal to e^^. Site 2 may be black (/i2 > h or 
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equivalent ly p2 < p) with an a priori probability p, or white with an a priori 
probability (1-p). If site 2 is white (/i2 < /i) , it must be down because white 
sites are yet to be relaxed. Thus, the conditional probability that the spin 
at site 2 is down, and the quenched field at site 2 is larger than h is given by 

Proh{2 i |1 lp2 <p) = [e-^ - {I-P)] ■ (16) 

The probability that the quenched field /i2 lies in the range h and h-dh, 
or equivalently p2 is in the range p and p+dp can be obtained by taking the 
derivative of the above expression. We obtain, 

Prob[p < p2 < p + dp]dp = [1 — e^^]dp (17) 

Similarly, 

Prob[p < pi < p + dp]dp = [1 — e^^]dp (18) 

We now address an issue which is crucial for determining ramp-II correctly. 
This concerns two adjacent doublets on plateau-I as shown in Figure (4). 
Let /ii, /i2, /is, h^ and h^ denote the quenched fields at sites 1,2, 3, 4, and 5 
respectively. If /i2 > /ii, and h^ > h^, then spins at sites 2,3, and 4 will be up 
at some value of the applied field on ramp-II. When this happens, i.e when 
a triplet of up spins is created on ramp-II, the central spin at site 3 becomes 
unstable and it fiips down. It stays down till the system reaches ramp-III. In 
order to take this effect into account, we must know the probability per site 
of observing two adjacent doublets on plateau-I, and also the distribution of 
fields at sites 2 and 4. 

12 



A doublet on plateau-I has an important property. It separates the lattice 
into two parts (one on each side of the doublet) which have evolved unin- 
fluenced by each other on ramp-I. Thus, we can separate Figure (4) into 
three parts as enclosed in the dashed boxes. Evolution inside each box has 
remained shielded from the outside on ramp-I. The evolution in the middle 
box requires that site 3 flips up before site 2 or site 4. The probability for 
this event is equal to |. Given this event, the probability that spins at sites 1 
and 5 remain down all along ramp-I is each equal to -. Thus the probability 
per site of observing two adjacent doublets on plateau-I is equal to ^. Note 
that it is quite different from the square of the probability of finding a single 
doublet! 

The shielding property of the dashed boxes in Figure (4) can also be 
used to calculate the a posteriori distribution of fields /ii, /i2, • • • , /is- The 
distribution of hi and h^ is the same as obtained above for a doublet with 
similar evolution on the two sides (Figure 3). 

Prob{p < pi < p + dp)dp = [1 — e^^](ij9 (19) 

and, 

Prob{p < J95 < p + dp)dp = [1 — e~'^]dp (20) 

We now turn to the distributions of /i2, and /i4. Suppose, in the middle 
box in Figure (4), the central spin at site 3 is up when all sites in the system 
with quenched fields hi > h (black sites) have been relaxed. This necessarily 
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means that h^ > h, i.e. site 3 is a black site, but sites 2 and 4 have other 
options. Site 3 will be up with probability |j9^ if 2 and 4 are both black, 
probability p{l — p)'^ if 2 and 4 are both white, and (1 — p)p'^ if 2 and 4 are 
mixed. Thus the probability of observing two adjacent doublets with the up 
spin separating them having field h^ > h (j93 < p) is given by 

1 



Piimip-s <p) 



-p^+p{i-pf + {i-p)p^ 



(21) 



In order to calculate the distribution of h^ and /i4, it is convenient to write 
the distributions of the smaller and the larger of these two fields separately. 
Without any loss of generality, we can assume 

/i2 = mm(/i2, h\) and h^ = vnaxiji^^ h\) 

If /i2 > /i, then we have h^ > h^ > h as well. Thus the fraction of adjacent 
doublets on plateau-I with /i2, /is, and h^ all greater than an arbitrary value 
h is given by 



p^ 1 



PiiUiiPs < Pi < P2 < p) = ^:^ (22) 



The above equation gives the cummulative fraction of j92 ^ P sites. The 
fraction of sites in the range p+dp > P2 ^ P can be obtained by differentiating 
the above expression. We get, 

PuniiP^P2<P + dp)=p^^ (23) 

The distribution of /i4 is obtained similarly. We find, 

niTii(l5 <PA<P + dp) = 2p{l -p)- (24) 
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Given that we have a pair of adjacent doublets, the probability densities 
of the quantities j9i, j92, Pa, and p^ (each normalised to unity) are given by: 



Pi{Pi) = e\l- e 



-pi 



(25) 



P2{P2) = ^pI 



(26) 



^4(^4) = 6^34(1 -;p4) 



(27) 



PbiPb) = e 



1 



-P5 



(28) 



In the following, we focus on adjacent doublets which create up triplets on 
ramp-II. These are the objects with h2 > hi and h4> h^. At this stage, we 
can determine the lower and the upper bound on these objects at any point 
on ramp-II when spins with hi > h have been relaxed. Suppose we order the 
adjacent doublets in the order of increasing /i2 or increasing h^. Note that a 
sequence in increasing /i2 does not posses any particular order in /i4, and vice 
versa. The lower bound is given by the fraction of objects with h > h2 > hi 
{h/i being free to have any value in the range h2 < h^ < A). We obtain 



p fi 

P2{p2)dp2 / pi{pi)dpi 

3 

= (6 + 6j9 + 3/)e^"^ + (1 + e)j9^ - -e/ - 6e 

23 1 
16 — — e (at j9 = 1) = .369 (approximately). 



(29) 



Thus a minimum of approximately 37 % adjacent doublets will give rise to 
(unstable) up triplets on ramp-II. 
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Similarly, the fraction of adjacent doublets with h > h4 > h^ > —A is 
given by 

/ p4{p4)dp4 f p5{P5)dp5 

3 

= (1 + e)/(3 - 2p) - 2e/ + -e/ + 6e - 6(1 + j9 + /)e^"P 



He - 17 
2 



^at j9 = 1) = .668 (approximately). (30) 



This gives the upper bound. No more than approximately 67 % of the adja- 
cent doublets can create (unstable) up triplets on ramp-II. 

VI. RAMP-II 

Ramp-II is determined by the combination of two opposite terms. The 
dominant term is the increase in magnetization due to the decrease in the 
number of doublets. When a doublet disappears, it adds an extra up spin in 
the system which increases the magnetization. Occasionally, a disappearing 
doublet creates a string of three up spins. A triplet of up spins is unstable on 
ramp-II, and the central spin of the triplet flips down as soon as the triplet 
is created. This decreases the magnetization. In the following, we calculate 
the above two terms separately. 

Refer to Figure (3) for calculating the first term. The probability that 
the doublet shown in the Figure (3) disappears when spins with hi > h are 
relaxed is given by 



^n = ^ i) p(P^)dPi J p{P2)dp2 (31) 
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The factor -\ is the probabihty per site of finding a doublet shown in Figure 
(3). The factor 2 takes care of the fact that either hi or h2 may flip up first. 
The integrals are written on the assumption that the spin at site 1 flips up 
first. Together they give the probabihty that hi > h, and h2 < hi. Note that 
when a doublet disappears, a pair of adjacent up spins is created. This is the 
reason for the choice of the subscript on P^Y- The superscript indicates that 
the probability refers to ramp-II. We obtain, 

P" = y [{^ + ^-') - {P + ^-'f (32) 

We now caculate the fraction of (unstable) up triplets formed in the system 
when all spins with hi > h have been relaxed on ramp-II. Let us refer to 
Figure (4). Recall that /i2 < /i4 in this figure. We want the probabihty that 
/i2 > /ii, /i4 > /i5, and /i2 > h. This is given by. 



Ptn = ^ Jo P^P^^'^P^ Lp(p^'>'^p^ 



Jo' p{PA)dpA 4 p{P5)dp5 



(33) 



/(?' p{PA)dp4 

The first factor is the probability per site of finding the object shown in Figure 
(4). The next two integrals give the probability that /i2 > /i, and hi < h2- 
The quantity in the square brackets is understood as follows: When /i2 is in 
the range h and h+dh, /i4 can be anywhere in the range /i2 to A. Let p{h/i) 
be the density of h/i in this range. Clearly, 

P2(fe2)=X^w(A4)§ (34) 



or. 
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Thus, 



or, 



P2{P2) = J^^ MP^dpA 



(35) 



Pa{P2) = 



dp2{p2) 

dp2 



(36) 



PtiPi) = 6P4 



(37) 



We get, 






3 + ^_6(l + i)p + 3p^+(l + i)p' 



ih^y^i^' 



- le (l + ^] - 6p - 3 (l + ^] p^ + 2pA e-P +(^ + 3p] e-^P 

We show in Figure (5) a comparison of the above expression with a result from 
the simulation. As may be expected, the agreement between the simulation 
and the theory is excellent. 

Incidentally, an interesting quantity is the ratio of (unstable) up triplets at 
the end of ramp-II to the fraction of adjacent doublets at the start of ramp- 
II. We noted in the previous section that this ratio must lie approximately 
in the range .369 to .668. If there were no correlations between adjacent 
doublets, this ratio would be equal to |, because the events /i2 > /ii, or 
/i4 > /i5 would occur with probability |. The exact value of the probability 
that /i4 > /is > —A and h2> hi> —A is given by 



^^^^ - 3e2 



11 7 



27 



\ — e e 

2 4 20 



(38) 



The quantity in the square brackets is approximately equal to .281. We have 
checked this result rather carefully numerically, and it is born out by the 
simulations. Finally, putting the various terms together, the probability of 
an up spin on ramp-II is given by 

P^iP) = ^/(i) + ^tY(^) - ^Tn(^) (39) 

The magnetization on ramp-II is given by 

m^\p) = 2PWp) - 1 (40) 

This expression has been superposed on the numerical data for ramp-II 
shown in Figure (1). The agreement between the numerical data and the 
theory is extremely good. The exact value of the magnetization on plateau- 
II, and its numerical estimate are given by 

"27 7 



m'\l) = 



30 6 3 



= .109 (approximately) (41) 



VII. PLATEAU-II 

Each down spin on plateau-II is a singlet. However, there are three dif- 
ferent classes of singlets: the singlets formed on ramp-I; singlets formed on 
ramp-II by a vanishing doublet; and finally the singlets formed on ramp-II 
by the unstable central spin of an up triplet flipping down. Each class is 
characterized by its own a posteriori distribution of the random field. 
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Let us denote the three distribution densities by p{^ , /o^^, and pl^ respec- 
tively. Let P/^, P2^ , and P^^^ denote the cummulative populations in each 
class when spins with hi > h have been relaxed on ramp-IIL We have, 

PF(p) = I' pl'{P^)dp^, i=l,2,3. (42) 

It is useful to think of the singlets in each class as being black (hi > h), or 
white {hi < /i), where hi is the quenched random field at the singlet site, and 
h is an arbitrary reference field. In order to calculate ramp-III, we need only 
the populations of black singlets in each class given by Pl^{p). If needed, 
one can obtain the density pl^{p) by differentiating Pl\p) with respect to p. 

The fraction of black singlets created on ramp-I is given by. 



1 



K b) = 1^ - 2 [1 - e^'1 - - [^'' - (1 - ^)] (43) 

The explanation of the above equation is as follows. Imagine ordering the 
sites of the lattice in order of decreasing quenched field on the site. When 
all sites with hi > h have been relaxed, the fraction of the relaxed sites is 
equal to p (the black sites). This fraction is made of the up sites (the second 
term on the right), black doublet sites (the last term), and the black singlets. 
Hence the equation for P(^(p). The last term is written as follows. In each 
doublet, there are two sites from which we can choose one. This accounts for 
the factor 2. The quantity in the square bracket gives the probability that 
the chosen site is black, and - is the probability that the other site can have 
any allowed value of the quenched field. 
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The fraction of black singlets generated by vanishing doublets on ramp-II 
is given by, 



j/// X r _„ /. Ni2 



Pi\p) = [e-' - (l - p)\ (44) 

The above equation is easily understood. It is the probability that both sites 
of the doublet are black. If both sites of the doublet are black, the one with 
higher random field must flip up on ramp-II, leaving us with a singlet on 
plateau-II which is black. 

The fraction of black singlets created by unstable triplets requires the 
calculation of triplets. We have calculated the fraction of triplets as they are 
formed on ramp-II. What we need now is a similar but different calculation. 
The point can be understood with a reference to Figure (4). Recall that h-^ > 
h/i > h^. On ramp-II, we needed the fraction of triplets with h^ > h, because 
the formation of triplets is controlled by this threshold. The restoration of 
the triplets on ramp-III is controlled by the condition h^ > h. In the earlier 
calculation, sites 2, 3, and 4 were all black sites. In the calculation needed 
now, only site 3 is black. Sites 2, and 4 are white, and we want hi < /i2, and 
^5 < ^4- The probability for this event is given by 



Pf{p) = 2 r dps /' dp, I' [1 - e-P'] dp, [' dp2 [' [1 - e-'P'] dp 



/O Jp3 Jpi 

As a check, we note that 



2 / dpA dp2 = (1 -p^y 



JP3 Jp4 
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Thus, the double integral gives the probability that sites 2 and 4 are white 
( compare with the second term on the right hand side of Equation 21 ). 
The extra terms in P^^ip) take into account the requirements h^ < h^, and 
hi < /i2. 



pfip) = 2jyp,i] 



l + e 



-1 



Pa 



-Pi 



dp4 



Pi 



l + e 



-1 



P2 



-P2 



dp2 



We get, 



pr(p) = 



1 + ^ + 



1.2, A 
4 ^ e ^ e2 



P 



1 + i 



/ + ip^ + 



1 + 



1 h + 5 , A] 2 , 1 [3 , 4 , X 



2^^- 



j3^ 



j3e-P + p'^e-P + I fl - 6-2^1 (45) 



As a check we note that. 
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1 
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This is the same as given by Equation (38). As one may expect, the number 
of triplets with h^ > h is larger than triplets with /i2 > /i, but the two merge 
dut h = —A i.e. p = 1. 

It is of some interest to calculate the number of triplets with h^ > h 
as well, although this quantity is not directly needed in the calculation of 
ramp-III. Denoting this quantity by P/^(j9), we get 



Pl\p) = 2 pdPA 



1 



-P5 



1 



dp5 / dp2 

•'Pa ''P2 



1 r 



1 



-Pl 



dpi 



Note that either h2 or h^ could have been the larger of the two fields 
although we have assumed that /i4 is the larger one. This accounts for the 
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factor 2 on the right hand side. The factor p4 comes from integrating p^ over 
the range to ^4. The other terms are self-explanatory. We get, 



PV{p) = 



1 + 21 , 1 N , 5 , A' 



P^ 



2 [3 , 4 , ^1 3 , 3 N , 11 4 



-Ip^+ [{l + l){l+p + p^)+p 



4 

i+P 



-2p 



(46) 



It can be checked that Pl^{l) also reduces to the expression in Equation (38) 
as may be expected. In Figure (5), we have shown the theoretical expressions 
for P^^{p), P^^{p), and Pl\p) along with the results from the simulations 
for the same quantities. The simulations were performed for a system of 
10^ spins, and averaged over 5-10^ different realizations of the random field 
distribution. The agreement is excellent as may be expected from an exact 
analytic result. The agreement between the simulation and the theory also 
justifies (albeit post facto) the implicit assumption in our analysis that the 
system is self-averaging. The fact that simulations over a relatively small 
size of the system (10^ spins) agree with the exact result is due to the super- 
exponential decay of correlations in this system 



VIII. RAMP-III 

We are not in a position to write the analytic expression for ramp-III 
at this stage. In the previous section, we have obtained some ingredients 
which are necessary for this purpose, but these are not sufficient. If the only 
process occurring on ramp-III were the gradual decrease in the number of 
singlets initially present on plateau-II, then we have the required information 
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to determine ramp-III. However, additional singlets are created (and later 
destroyed) on ramp-III. Let us illustrate this with an example. Consider a 
down spin on plateau-II, say with random field hd- Its nearest neighbors are 
necessarily up because plateau-II contains only singlets. The next nearest 
neighbors can be up or down. Consider a nearest neighbor, and let hu be 
the random field on it. The site with field hu may have one or both of its 
nearest neighbors down. Suppose one of the nearest neighbors of hu is up. 
Now consider the situation when the singlet just fiips up on ramp-III at an 
applied field ha, i.e. —2\J\ + hd + ha = e, where e > 0. The applied field at 
this point is /ia = 2| J| — /i^ + e. If hu is to remain up after hd has fiipped up, 
we must have hu — 2\J\ + ha> 0, or hu — hd + e > 0. This is not guaranteed 
if hu < hd- In this case hu will fiip down when hd fiips up, and later fiip 
up again in a higher applied field. The conditional probability distributions 
which govern the creation of these new singlets and their restoration on ramp- 
III require further careful investigation. We hope to take up this work in a 
future article. 

IX. CONCLUSION 

We have presented an exact solution of the zero-temperature dynamics of 
the one-dimensional anti-ferromagnetic random field Ising model in a slowly 
increasing applied field for a limited range of the applied field. This problem 
was posed a few years ago, but to the best of our knowledge remained un- 
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solved in this period. The solution presented here complements the analytic 
solution of the non-equilibrium dynamics of the ferromagnetic random field 
Ising model. The key difference from the ferromagnetic case is the absence 
of avalanches. A large number of physical systems relax by avalanches, but 
clearly not all of them. Therefore, models of relaxation without avalanches 
have to be studied as well. The analysis presented here is a small step in this 
direction. It also attempts to develop the techniques used in the study of the 
problem of random sequential adsorption (RSA). The relationship between 
ramp-I and RSA was noted earlier. Indeed, the RSA rate equations were 
used to solve the problem of ramp-I. However, the problem of ramp-II and 
ramp-III could not be solved. Here we have solved the problem of ramp-I 
by an alternate (although not profoundly different) method which focuses 
on the a posteriori modification of the random field distribution by the dy- 
namical selection process. This method yields more detailed information on 
ramp-I than was available earlier, and allows us to get an exact expression 
for ramp-II. We hope that the method employed here will enable us to solve 
the problem of ramp-III in the near future, and it may be useful in other 
related problems as well. 

We thank D Dhar for useful discussions during a visit to NEHU. 
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Figure Captions 

Figure 1: 
Magnetization m{ha) in an applied field ha for an anti-ferromagnetic 
RFIM (J = —1), obtained from a simulation of 10^ spins aver- 
aged over 10^ different realizations of the random field distribution 
of width 2 A = 1. The solid line shows the magnetization in increas- 
ing field. The broken line shows the return half of the hysteresis 
loop. The data in increasing field is separated into five parts along 
the applied field: Ramp-I {ha = — 2| J| — A to ha = — 2| J| + A); 
Plateau-I {ha = -2| J| + A to /la = -A); Ramp-II {ha = -A to 
ha = +A); Plateau-II {ha = +A to /la = 2| J| — A); and Ramp-III 
(/la = 2| J| — A to /la = 2| J| + A). Theoretical expressions obtained 
in the paper have been superposed on Ramp-I, Plateau-I, Ramp-II, 
and Plateau-II in increasing field. 

Figure 2: 
Spins on Ramp-I in an appfied field — 2| J| — h. Filled circles show 
sites with quenched field hi > h. The probability per site of a doublet 
(two adjacent down spins) such as AB is equal to e~^^, where p is 
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the fraction of filled circles on the infinite lattice. 

Figure 3: 
A doublet on Plateau-I: hi and /12 are the quenched random fields 
on the doublet sites 1 and 2 respectively. 

Figure 4: 
Two adjacent doublets on Plateau-I: Each doublet separates the lat- 
tice into two parts whose evolution histories on Ramp-I are indepen- 
dent of each other. Evolutions inside each dashed box is shielded 
from outside. The probability that spin at site 3 fiips up on Ramp-I 
is therefore equal to |. Given this, the probability that the spins 
at sites 1 and 5 remain down all along Ramp-I is equal to - each. 
The shielding property of the boxes can also be used to determine a 
posteriori distribution of random fields /ii, /125 ^3, ^4, and /15. 

Figure 5: 
Unstable up triplets on Ramp-II with h2 > h (lower curve); h4> h 
(middle curve); and hs > h (upper curve). Refer to Figure 4 for /12, 
/14, and /13. The simulation data is shown by lines, and the symbols 
on each line show the corresponding theoretical prediction. 
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